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Introduction 


Single crystal superalloy turbine blades used in high pressure turbomachinery are subject to 
conditions of high temperature, triaxial steady and fatigue stresses, fretting stresses in the blade 
attachment and damper contact locations, and exposure to high-pressure hydrogen. The blades 
are also subjected to extreme variations in temperature during start-up and shutdown transients. 
The most prevalent HCF failure modes observed in these blades during operation include 
crystallographic crack initiation/propagation on octahedral planes, and noncrystallographic 
initiation with crystallographic growth. Numerous cases of crack initiation and crack propagation 
at the blade leading edge tip, blade attachment regions, and damper contact locations have been 
documented. Understanding crack initiation/propagation under mixed-mode loading conditions 
is critical for establishing a systematic procedure for evaluating HCF life of single crystal turbine 
blades. 

Techniques for evaluating two and three dimensional subsurface stress fields in anisotropic 
contacts are presented in this report. Figure 1 shows typical damper contact locations in a turbine 
blade. The subsurface stress results are used for evaluating contact fatigue life at damper contacts 
and dovetail attachment regions in single crystal nickel-base superalloy turbine blades. 



Damper 

Contact 


Figure 1 : Damper contact locations on a typical turbine blade 
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Two-Dimensional Stress Distribution in an Anisotropic Elastic Half-Space 


The damper contact regions shown in Figure 1 will be modeled as an elastic anisotropic half- 
space. This approximation is reasonable since Hertzian type contact stresses are confined to very 
small volumes in the vicinity of the contact. An analytical procedure will be presented for 
evaluating the subsurface stresses in the elastic half-space using a complex potential method 
outlined by Lekhnitskii [1]. Figure 2 shows the elastic half-space subjected to normal traction 
N(£) and tangential traction T (^ ) over the region -a to +a on the x-axis. The traction forces are 
independent of z, and functions of x and y only. The stresses are also functions of x and y only. 
Under the assumptions of generalized plane strain, the subsurface stresses due to the applied 
traction forces can be determined as outlined below. 

The stress functions are given by 
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The jUi are the roots of the cylindrical characteristic equation, given by Eq. (4), and z = x + fly. 
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Figure 2: Anisotropic elastic half-space under generalized plane deformation 
subjected to normal and tangential traction forces. 
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The matrix a l} relates the strains to the stresses. The a \ j are functions of the crystal orientation. 
The stresses are then given by: 
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The normal traction force N(£) is given by the Hertzian cylindrical contact force as 
= Po ^~e/a 2 and T (g) = ju f p 0 ^l-£ 2 /a 2 , where p 0 is the peak pressure and // f the 


coefficient of friction. The stress solution has been programmed in Mathcad and subsurface 
stresses computed for various crystal orientations. Figure 3 shows a representative half-space <7 y 
stress distribution for a = 0.01 inch, po = 260 ksi, and for the (x, y, z) axes parallel to the edges 
of the FCC crystal, i.e. x = <100>, y = <0 1 0> and z = <001>. The analytical solution shows 
good agreement with the finite element numerical solution. 


Subsurface fretting fatigue cracks at the damper contacts caused by oscillatory tangential 
stresses, T(^), is of much practical interest. Figure 4 shows the variation in resolved shear stress 

Finite Element Solution 



Analytical Solution 



Figure 3: Stress (<7 y ) contours using analytical solution and finite element (Ansys) solution 
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amplitude on the primary octahedral planes, as a function of secondary crystallographic 
orientation. It is seen that the highest resolved shear stress amplitude, Al rss , varies 32 percent as 
the secondary orientation varies from 0 to 90 deg. The secondary orientation of the blades is not 
controlled during the casting process and hence its variation can result in blade -to-blade variation 
in fatigue stresses (Al rss ) at constant load, which is clearly undesirable from a design standpoint. 

Conclusions 


The two-dimensional analytical solution presented is an efficient and accurate method for 
obtaining contact stresses in anisotropic half-spaces. Finite element solution of contact problems 
is very tedious and requires highly refined meshes to obtain accurate stress solutions. Variation 
in secondary orientation between blades can result in large variation in fatigue life at damper 
contact locations. 

Three-dimensional contact stress solutions for anisotropic half-spaces has also been obtained for 
elliptical contacts. Stress solutions have been obtained by linear superposition of a point-load 
solution outlined in Lekhnitskii [1]. The details of this solution could not be included in this 
report because of lack of space. Complete details of the 2D and 3D solutions will be documented 
in a future NASA Technical Publication. 
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Figure 4: Effect of secondary orientation on resolved shear stress amplitude (Al rss ) 
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